Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 February 2008 (MN IATpX style file v2.2) 



The SAURON project - IX. A kinematic classification for early-type 
galaxies. 

Eric Emsellem, 1 Michele Cappellari, 2 ' 3 Davor Krajnovic, 3 Glenn van de Ven, 2 ' 4 ' 5 * R. Bacon, 
M. Bureau, 3 Roger L. Davies, 3 R T. de Zeeuw, 2 Jesus Falcon-Barroso, 2,6 
Harald Kuntschner, 7 Richard McDermid, 2 Reynier F. Peletier, 8 Marc Sarzi, 9 

1 Universite de Lyon, France; universite Lyon 1, F-69007; CRAL, Observatoire de Lyon, F-69230 Saint Genis Laval; CNRS, UMR 5574 ; ENS de Lyon, France 
2 Sterrewacht Leiden, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands 

3 'Sub-Department of Astrophysics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, United Kingdom 
^Department of Astrophysical Sciences, Peyton Hall, Princeton, NJ 08544, USA 
5 Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 

6 European Space and Technology Centre (ESTEC), Keplerlaan 1, Postbus 299, 2200 AG Noordwijk, The Netherlands 

7 Space Telescope European Coordinating Facility, European Southern Observatory, Karl-Schwarzschild-Str 2, 85748 Garching, Germany 

^Kapteyn Astronomical Institute, Postbus 800, 9700 AV Groningen, The Netherlands 

® Centre for Astrophysics Research, University of Hertfordshire, Hatfield, Herts ALIO 9AB 



4 February 2008 



ABSTRACT 

Two-dimensional stellar kinematics of 48 representative E and SO galaxies obtained with 
the SAURON integral-field spectrograph reveal that early-type galaxies appear in two broad 
flavours, depending on whether they exhibit clear large-scale rotation or not. We define a new 
parameter \r = (R \V\) / (R \/V 2 + a 2 }, which involves luminosity weighted averages over 
the full two-dimensional kinematic field, as a proxy to quantify the observed projected stellar 
angular momentum per unit mass. We use it as a basis for a new kinematic classification: 
early-type galaxies are separated into slow and fast rotators, depending on whether they have 
Afl values within their effective radius R e below or above 0.1, respectively. Slow and fast ro- 
tators are shown to be physically distinct classes of galaxies, a result which cannot simply be 
the consequence of a biased viewing angle. Fast rotators tend to be relatively low luminosity 
galaxies with Mb > —20.5. Slow rotators tend to be brighter and more massive galaxies, 
but are still spread over a wide range of absolute magnitude. Three slow rotators of our sam- 
ple, among the most massive ones, are consistent with zero rotation. Remarkably, all other 
slow rotators (besides the atypical case of NGC 4550) contain a large kpc-scale kinematically 
decoupled core (KDC). All fast rotators (except one galaxy with well-known irregular shells) 
show well aligned photometric and kinemetric axes, and small velocity twists, in contrast with 
most slow rotators which exhibit significant misalignments and velocity twists. These results 
are supported by a supplement of 1 8 additional early-type galaxies observed with SAURON. In 
a companion paper (Paper X), we also show that fast and slow rotators are distinct classes in 
terms of their orbital distribution. We suggest that gas is a key ingredient in the formation and 
evolution of fast rotators, and that the slowest rotators are the extreme evolutionary end point 
reached deep in gravitational potential wells where dissipationless mergers had a major role in 
the evolution, and for which most of the baryonic angular momentum was expelled outwards. 
Detailed numerical simulations in a cosmological context are required to understand how to 
form large-scale kinematically decoupled cores within slow rotators, and more generally to 
explain the distribution of A r values within early-type galaxies and the distinction between 
fast and slow rotators. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: forma- 
tion - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 



Hubble Fellow 



The origin of the classification fork for galaxies can be found in 
an early paper by Jeans (1929), with the SOs as a class being in- 
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troduced by ? to account for the important population of flattened 
objects in nearby clusters (?). In a recent debate on galaxy classi- 
fication, ? mentioned that the simplest definition of an SO galaxy 
remains "a disc galaxy more flattened than an E6 elliptical but 
with no trace of spiral arms or recent star formation". Elliptical (E) 
and lenticular (SO) galaxies are usually gathered into the so-called 
early-type category, and are recognised to share a number of global 
properties (?) such as their relatively low dust and interstellar gas 
content and their overall red colours. The Hubble sequence is how- 
ever seen as a continuous one between ellipticals and spirals, with 
the SOs occupying the transition region with typical bulge to disc 
ratios of ~ 0.6. SOs are thus considered disc-dominated galaxies, 
while Es are spheroid-dominated. 

Such contrived galaxy types may be misleading, most evi- 
dently because "the sequence E0-E6 is one of apparent flattening" 
(?, hereafter KB96). A modern classification scheme should go 
beyond a purely descriptive tool, and should therefore encompass 
part of our knowledge of the physical properties of these objectfl 
This was advocated by KB96 who wished to update the Hubble 
sequence by sorting ellipticals in terms of the importance of ro- 
tation for their stellar dynamical state. They used the disciness or 
boxiness of the isophotes to quantify anisotropy and to define re- 
fined types: E(d) galaxies (for discy ellipticals) making the link be- 
tween SOs and E(b) galaxies (for boxy ellipticals). The disciness 
(or boxiness) was then provided by a measure of the now classi- 
cal normalised a^/a term (see ?, for details): positive and negative 
cn/a terms correspond to discy and boxy deviations from ellipses, 
respectively. 

This extension of the Hubble types has the merit of upgrad- 
ing our view of Es and SOs via some easily accessible observable 
parameter, and it follows the philosophy that a mature classifica- 
tion scheme should include some physics into the sorting criteria. 
It does, however, use a photometric indicator as an attempt to quan- 
tify the dynamical state of the galaxy, which may be unreliable. 
More importantly, it conserves the dichotomy between SOs and Es, 
relying on the old (and ambiguous) definition of an SO. 

We have recently conducted a survey of 72 early-type (E, 
SO, Sa) galaxies using the integral-field spectrograph SAURON 
mounted on the William Herschel Telescope in La Palma (?, here- 
after Paper I; ?, hereafter Paper II). This allowed us to map the stel- 
lar and gas kinematics as well as a number of stellar absorption line 
indices up to about one effective radius R e for most of the galaxies 
in the sample. The two-dimensional stellar kinematics for the 48 
E and SO galaxies (?, hereafter Paper III) show a wide variety of 
features such as kinematically decoupled or counter-rotating cores, 
central discs and velocity twists. More importantly, there seem to be 
two broad classes of observed stellar velocity fields, with galaxies 
in one class exhibiting a clear large-scale rotation pattern and those 
in the other showing no significant rotation (Fig.[T](. The existence 
of these two classes must be linked to the formation and evolution 
of early-type galaxies, and is in any case a key to understand their 
dynamical state (see also ?, hereafter Paper IV). 

Using the unique data set obtained in the course of the 
SAURON project, we here revisit the early-type galaxy classi- 
fication issues mentioned above, using the available full two- 
dimensional kinematic information. A companion paper (?, here- 
after Paper X) examines in more detail the orbital anisotropy of el- 
liptical and lenticular galaxies using the same data set. After a brief 
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presentation of the data set and methods (Sect.[2]l, we define a new 
parameter as a proxy to robustly quantify the angular momentum 
of galaxies (Sect.[3](- In Section[4] we examine how this parameter 
relates both to more standard photometric classification schemes 
including de Vaucouleurs morphological classification of Es and 
SOs (?), or the revision of the Hubble sequence by KB96, and kine- 
metric properties of early-type galaxies (Sect. [4~4b . We then briefly 
discuss the implications of our results on the potential scenarios for 
the formation and evolution of these galaxies (Sect.O, and con- 
clude in Section[6] 



2 DATA AND METHODS 

2.1 The SAURON sample of E and SOs 

The SAURON sample has been designed as a representative sam- 
ple of 72 nearby (cz < 3000km s -1 ) early- type galaxies in the 
plane of ellipticity e versus absolute B band magnitude Mb- We 
restricted the sample to 24 objects for each of the E, SO and Sa 
classes, with 12 'cluster' and 12 'field' targets in each group. Al- 
though this sample is itself not a complete one, it has been drawn 
in a homogeneous way from a complete sample of more than 300 
galaxies (see Paper II for more details about the full SAURON sam- 
ple). We will therefore solely focus on the SAURON E+S0 sample 
(which thus contains 48 galaxies) for conducting our analysis in 
Sections[2]to|4] An additional set of 18 E+S0 galaxies has been in- 
dependently observed with a similar SAURON setup in the course of 
various other projects. These objects will be treated as "specials", 
most having features which motivated a specific observation, and 
will be only briefly mentioned in Sect. [5] to confront the obtained 
results. 



2.2 Photometry 

Ground-based photometric MDM Observatory data (Falcon- 
Barroso et al. in preparation) were obtained for all galaxies of the 
SAURON sample. We also made use of additional Hubble Space 
Telescope WFPC2 data which are available for 42 galaxies out of 
the 48 E/S0. SAURON reconstructed images were used as well to di- 
rectly derive a global ellipticity e via moments of the surface bright- 
ness images as: 



1 



<2/!> 
<* 2 > 



(1) 



where x and y are the sky coordinates along the photometric major 
and minor axes respectively (see Paper X for details). The brackets 
( ) correspond to a flux weighted sky average as in: 



<G> = 



(2) 



where the sums extends over a predefined region on the sky, with 
N p the number of pixels, and Fi the flux value within the i" 1 pixel. 
In the following Sections, we denote e and e 1( / 2 as the global ellip- 
ticity computed from the SAURON data within 1 R e and R e /2 or 
restricted to the equivalent effective aperture of the SAURON field 
of view, whichever is smaller. These will be used in Sect. [3] where 
we examine kinematic quantities derived from SAURON data as a 
function of the flux weighted global ellipticity. 

We also obtained radial profiles for the ellipticity and a^/a 
(with a the ellipse semi major-axis) parameters from the MDM and 
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HSTAVFPC2 photometric data sets using the GALPHOT pack- 
age (?). We then derived e e and (0,4, / a) e , the mean ellipticity and 
0,4 /a values within 1 R e , by taking the flux weighted average of 
the GALPHOT profiles taking into account the corresponding sky 
area. More specifically, the mean G e of a quantity G(R) within 
1 R e derived from its radial profile, where R is the semi major-axis 
radius, is defined as (see e.g. ?) 

G _ J*'q(R)F(R)G(R)RdR 
f*°q(R)F(R)R6R ' 

where q(R) and E(-R) are respectively the best-fit ellipse axis ratio 
and surface brightness profiles. Using the sampled radial profiles, 
we approximate this with 

Eti g(R k )F(R k )G(R k ) (R 2 ouUk - Rf n k ) 
£ ~ Ef=i q(Rk)F(R k ) (Rl ut k - Rf n k ) ' 

where 7?j n k and -Rout.fc correspond to the inner and outer radii 
of the k' h annulus. The values for e, e e and (o4/a) e are provided 
in Table [T] Observed differences between e and e e values are al- 
most always smaller than 0. 1 , and due to their respective luminos- 
ity weighting, e is measured directly from the SAURON data, and 
will therefore be confronted with other quantities measured from 
the same data set, while e e being derived from one-dimensional ra- 
dial profiles can serve in comparison with previous works. For both 
e e and (0.4 /a) e , the overall agreement between our values and pub- 
lished ones (?) is excellent (see Paper X for a detailed comparison). 

Results from fitting Sersic and Sersic-core laws (??) to the 
radial luminosity profiles will be presented in detail in a subsequent 
paper of this series (Falcon-Barroso et al., in preparation). In the 
present paper, we will only mention trends (Sect. POt , considering 
the Sersic index n (where n — 1 corresponds to an exponential 
luminosity profile, and n = 4 to a de Vaucouleurs R 1 ^ 4 law), as 
well as the classification of the central photometric profiles with 
either shallow or steep inner cusps, labelled respectively as "cores" 
and "power-laws" (????). 

2.3 The SAURON data 

SAURON is an integral-field spectrograph built at Lyon Observa- 
tory and mounted since February 1999 at the Cassegrain focus of 
the William Herschel Telescope. It is based on the TIGER concept 
(?), using a microlens array to sample the field of view. Details of 
the instrument can be found in Paper I and II. All 48 E and SO 
galaxies were observed with the low resolution mode of SAURON 
which covers a field of view of about 33"x4l"with OV 94x0'.' 94 
per square lens. Mosaicing was used to cover up to a radius of 1 
R e . Only for the two galaxies with the largest R e (NGC4486 and 
NGC 5846), do we reach a radius of ~ 7? e /3 only. 

All data reduction was performed using the dedicated 
XSauron software wrapped in a scripted pipeline (Paper II). For 
each target, individual datacubes were merged and analysed as de- 
scribed in Paper III, ensuring a minimum signal-to-noise ratio of 60 
per pixel using the binning scheme developed by ?. The SAURON 
stellar kinematics were derived using a penalised pixel fitting rou- 
tine (?), which provides parametric estimates of the line-of-sight 
velocity distribution (hereafter LOSVD) for each spaxel. In Pa- 
per III, we have presented the corresponding maps, which include 
the mean velocity V, the velocity dispersion a and the Gauss- 
Hermite moments /13 and /14, forthe 48 E and SO SAURON galaxies. 

As mentioned in Paper III, these quantities were measured fit- 
ting all V, a, /13 and /14 simultaneously: this ensures an optimal 



representation of the corresponding LOSVD. In the present paper, 
we focus on the first two true velocity moments, ^1 and pi, which 
are sometimes estimated by use of the Gauss-Hermite expansion 
of the LOSVD (?). As emphasised in Paper IV, the second order 
velocity moment is very sensitive to the details of the high veloc- 
ity wings, which can rarely be accurately measured. We therefore 
decided to rely on a simpler but more robust single gaussian fit 
(excluding higher order velocity moments), and used the gaussian 
mean V and standard deviation a to approximate the first and sec- 
ond velocity moments. 

2.4 Kinemetry 

Following ? who recently advocated the use of a method gener- 
alising the isophotal-shape tools (???), we employ kinemetry as 
a quantitative approach to analyse the SAURON stellar kinematic 
maps. Applying kinemetry on a velocity map provides radial pro- 
files for the kinematic position angle PAkin, axis ratio q/cm, and 
Fourier kinemetry terms, the dominant term ki representing the 
velocity amplitude. We define a kinematic component to have con- 
stant or slowly varying PAkin and qtin radial profiles (taking into 
account the derived error bars). We then identify two separate com- 
ponents when we observe an abrupt change either with Ag^ > 0.1, 
or APAfci n > 10°, or a double-hump in fci with a local minimum 
in between. The transition between the two radial ranges is often 
emphasised by a peak in the higher order fcs Fourier term, which 
thus serves as an additional signature for such a change (?). Veloc- 
ity maps which exhibit the presence of at least two stellar velocity 
components are tagged as Multiple Component (MC), as opposed 
to Single Component (SC). 

We use a number of terms to represent some basic proper- 
ties of the individual kinematic components via quantitative crite- 
ria, following the definitions provided in ? : 

• Low-level velocity (LV): defined when the maximum velocity 
amplitude ki is lower than 15 km s _1 . Note that when the velocity 
amplitude is constant over the field, "PAkin and qtin are ill-defined. 
When the central kinematical component is LV, we label the galaxy 
as a central low-level velocity (CLV) system. 

• Kinematic misalignment (KM): defined when the absolute dif- 
ference between the photometric and kinemetric position angles 
(PAp ho t and PA kin ) is larger than 10°. 

• Kinematic twist (KT): defined by a smooth variation of the 
kinematic position angle PAfci„ with an amplitude of at least 10° 
within the extent of the kinematic component. 

A kinematically decoupled component (KDC) is defined as an 
MC having either an abrupt change in PAkin, with a difference 
larger than 20° between two adjacent components, or an outer LV 
kinematic component (which prevents that component to have a ro- 
bust PAkin measurement). This definition roughly corresponds to 
the more standard appellation of KDC used in the past (?), although 
the two-dimensional coverage provided by integral-field spectro- 
graphs allows a more sensitive detection procedure. These kine- 
metric groups will be examined in Sect. 14.41 



3 QUANTIFYING THE ANGULAR MOMENTUM 

The velocity fields presented in Paper III revealed a wealth of 
structures such as decoupled cores, velocity twists, misalignments, 
cylindrical or disc-like rotation (see also Section |4~4l l. It is difficult 
to disentangle the relative contributions of a true variation in the 
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Figure 1. SAURON Stellar velocity fields for our 48 E and SO galaxies (see Paper III), the global outer photometric axis being horizontal. Colour cuts were 
tuned for each individual galaxy as to properly emphasise the observed velocity structures. A representative isophote is overplotted in each thumbnail as a 
black solid line, and the centre is marked with a cross. Galaxies are ordered by increasing value of A^j c (from left to right, top to bottom; see Sect. [3). Slow 
rotators are galaxies on the first two rows. NGC numbers and Hubble types are provided in the lower-right and upper-right corners of each panel, respectively. 
Tickmarks correspond to 10 /; . 
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internal dynamical state and the effect of projection. A first exami- 
nation of the stellar velocity maps for the 48 SAURON E/SO galaxies 
(Fig. H] where velocity cuts have been adjusted as to properly em- 
phasise the observed velocity structures) suggests that early-type 
galaxies come in two broad flavours: one which exhibits a clear 
large-scale and rather regular rotation pattern, and another which 
shows perturbed velocity structures (e.g. strong velocity twists) or 
central kinematically decoupled components with little rotation in 
the outer regions. One way to constrain the internal dynamics for 
a specific galaxy is to build a detailed model using all available 
observables. This was successfully achieved in Paper IV on a sub- 
sample of 24 galaxies, for which accurate distances and high spatial 
resolution photometry are available, and where no strong signature 
of non-axisymmetry is observed. From these models, a consider- 
able amount of detail on the orbital structure was derived (see also 
Paper X). We, however, also need a simple measurable parameter 
which quantifies the (apparent) global dynamical state of a galaxy, 
and which is applicable to all galaxies in our representative sample 
of Es and SOs: this is the purpose of this section. 

3.1 From V/a to Xr : a new kinematic parameter 

The relation of V/a versus mean ellipticity e (hereafter the 
anisotropy diagram) was often used in the past to confront the ap- 
parent flattening with the observed amount of rotation (???). In the 
now classical treatment of the anisotropy diagram, the maximum 
observed rotational velocity V max and the central dispersion ao are 
generally used as surrogates for the mass-weighted mean of the 
square rotation speed and the random velocity. This was mostly 
constrained by the fact that stellar kinematics were available along 
at most a few axes via long-slit spectroscopy. ? has recently revis- 
ited this diagram to design a more robust diagnostic of the veloc- 
ity anisotropy in galaxies using two-dimensional kinematic infor- 
mation. Starting from the Tensor Virial Theorem, Binney reformu- 
lated the ratio of ordered versus random motions in terms of inte- 
grated quantities observable with integral-field spectrographs such 
as SAURON, namely (V 2 ) and {a 2 }, where V and a denote re- 
spectively the observed stellar velocity and velocity dispersion, and 
the brackets correspond to a sky averaging weighted by the surface 
brightness (see Eq.[2}. 

The SAURON data provide us with a unique opportunity to de- 
rive for the first time a robust measurement of V/a for a sample of 
local early-type galaxies. We have therefore derived (V 2 ) and (a 2 ) 
up to ~ lR e , the resulting V/a values for the 48 SAURON E and 
SO galaxies being given in Table[T] SOs and Es tend to have on av- 
erage relatively high and low V/a, respectively. A very significant 
overlap of the two types still exists, SOs reaching values of V/a as 
low as 0.1, and Es as high as 0.7: as expected, the Hubble type can 
obviously not serve as a proxy for the galaxy kinematics (see the 
distribution of Hubble types in Fig.[T]l. 

The availability of two-dimensional stellar kinematics makes 
V/a a useful tool to examine the dynamical status of early-type 
galaxies (see Paper X). It fails, however, to provide us with a way 
to differentiate mean stellar velocity structures as different as those 
of NGC 3379 and NGC 5813 (see Fig.Q). These two galaxies both 
have V/a ~ 0.14, and ellipticities e in the same range (0.08 and 
0.15, respectively), but their stellar velocity fields are qualitatively 
and quantitatively very different: NGC 3379 displays a regular and 
large-scale rotation pattern with a maximum amplitude of about 
60 km s , whereas NGC 5813 exhibits clear central KDC with a 
peak velocity amplitude of ~ 85 km s _1 and a mean stellar veloc- 
ity consistent with zero outside a radius of ~ 12". The derivation 



of V / a includes a luminosity weighting that amplifies the presence 
of the KDC in NGC 5813. As a consequence, NGC 3379 with its 
global rotation and NGC 5813 with its spatially confined non-zero 
velocities end up with a similar V/a. 

We therefore need to design a new practical way to quantify 
the global velocity structure of galaxies using the two-dimensional 
spatial information provided by integral-field units. The ideal tool 
would be a physical parameter which captures the spatial informa- 
tion included in the kinematic maps. Since we wish to assess the 
level of rotation in galaxies, this parameter should follow the na- 
ture of the classic V/a: ordered versus random motion. A measure 
of the averaged angular momentum L — (R A V), could play the 
role of V, and should be able to discriminate between large-scale 
rotation (NGC 3379) and little or no rotation (NGC 5813). Such a 
quantity, however, depends on the determination of the angular mo- 
mentum vector direction, which is not an easily measured quantity. 
We therefore use a more robust and measurable quantity {R \V\), 
where R is the observed distance to the galactic centre, as a surro- 
gate for L, and the brackets ( ) correspond to a luminosity weighted 
sky average (see Eq.[2}- We can then naturally define a dimension- 
less parameter, after normalising by e.g. mass, which leads to a 
proxy for the baryon projected specific angular momentum as: 



X f 



(R\V\) 



{RVV 2 +a 2 ) 
and measured via two-dimensional spectroscopy as: 



(5) 



(6) 



where Ft is the flux inside the i bin, Ri its distance to the cen- 
tre, and Vi and ai the corresponding mean stellar velocity and ve- 
locity dispersion. The normalisation by the second velocity mo- 
ment V 2 + a 2 implies that Xr goes to unity when the mean stellar 
rotation (V) dominates, with V 2 + a 2 being a reasonable proxy 
for mass (see Appendix [A}. The use of higher order moments of 
either V or the spatial weighting R would make this parameter 
more strongly dependent on the aperture and presence of noise in 
the data. Xr obviously depends on the spatial extent over which 
the sums in Eq. [6] are achieved. In practice, we measure XR(R m ) 
within regions defined by the photometric best fit ellipses, where 
R m is the mean radius of that ellipse (a VI — e, with a its semi- 
major axis and e its ellipticity), the area A e ai p3e of the correspond- 
ing aperture being thus tyR 2 ^, the area of a circle with a radius of 
R m - When our SAURON kinematic measurements do not fully sam- 
ple the defined ellipse, we instead set the radius R m = \J As/n 
as that of a circular aperture with the same area As on the sky 
actually covered by the SAURON data inside that aperture (see Pa- 
per IV). For a specific galaxy, we measure \R(R m ) up to the ra- 
dius for which we reach a maximum difference of 15% between 
As and A e iu pse : this guarantees that the SAURON kinematic data 
still properly fill up the elliptic aperture defined by the photometry. 



3.2 Rotators and A^ 

Values of A.r c , Xr within 1 R e or the largest radius allowed by our 
SAURON data, whichever is smaller, are provided in Table[TJfor the 
48 E/S0 galaxies. As expected, NGC 3379 and NGC 5813 which 
have similar V/a values (see Fig.[T]and Sect. |3. It but qualitatively 
different velocity structures, are now quantitatively distinguished 
with Xr c values of 0.14 and 0.06, respectively. This difference is 



6 E. Emsellem et al. 




R/Rc 



Figure 2. Radial Xr profiles for the 48 E and SO galaxies of the SAURON 
sample. Profiles of slow and fast rotators are coloured in red and blue, re- 
spectively. NGC numbers are indicated for all fast rotators and most slow 
rotators (a few were removed for legibility). 



significant because the formal uncertainty in the derivation of Xr c 
due to the presence of noise in the data is almost always negligible, 
and anyway smaller than 0.02 for galaxies such as NGC 3379 and 
NGC 5813 (Appendix [B}: this can be understood because Xr b in- 
cludes averages over a large area. There is, however, a systematic 
(positive) bias which obviously increases as the the velocity ampli- 
tudes in the galaxy decrease, and can reach up to about 0.03 in the 
measurement of Xr^ . This bias is therefore dominant for the three 
galaxies with very low \r c (< 0.03), the mean stellar velocities 
being in fact consistent with zero values everywhere in the field of 
view. 

As we go from galaxies with low to high Xr f values, the 
overall velocity amplitude naturally tends to increase. More impor- 
tantly, there seems to be a qualitative change in the observed stellar 
velocity structures. This is already illustrated in Fig.[T] where the 48 
SAURON stellar velocity fields are ordered, from left to right, top to 
bottom, by increasing value of Xr c . Rotators with Xr c < 0.1 ex- 
hibit low stellar mean velocities at large radii, with very perturbed 
stellar kinematics and large-scale kinematically decoupled compo- 
nents (this point will be further examined in Sect. 14.41 ). 

This qualitative change is nicely illustrated in Fig.|2]where we 
show the radial Xr profiles. Galaxies with Xr c below and above 
0.1 exhibit qualitatively very different Xr profiles: the former have 
either decreasing or nearly flat (and small amplitude) Xr profiles, 
while the latter preferentially exhibit significantly increasing Xr ra- 
dial profiles. Observed A_r gradients are therefore negative or rather 
small within 1 R e for galaxies with A_r c < 0.1, and we will label 
them as "slow rotators". This contrasts with the significantly rising 
A_r profiles for galaxies with Xr,, > 0.1, which all exhibit clear 
large-scale and relatively regular rotation patterns, and which we 
label, by opposition, as "fast rotators": the fact that this class in- 
cludes both mild and very fast rotators is discussed in Sect. |3.4| 

As mentioned, Table[T]includes Xr c values derived using the 
SAURON two-dimensional kinematic maps available and a default 
equivalent aperture of 1 R e . This aperture is in fact covered by the 
SAURON datacubes for 17 galaxies out of the 48 in the SAURON 
E/S0 sample (see Paper IV), with two galaxies being mapped only 
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Figure 3. Histogram of Xr: bars with plain colours correspond to Xr c = 
Xr(Rc) and the dashed bars to Xr(R c /2). Red bars indicate the Xr bin 
for slow rotators: the number of slow rotators does not change if we were 
to measure A^ at R c /2. 

to ~ 0.3 R e (NGC 4486 and NGC 5846). Galaxies with the nar- 
rowest relative spatial coverage (with the largest i? e : NGC 4486, 
NGC 5846) are among the slowest rotators of our sample: these 
two galaxies are in fact known not to exhibit any significant rota- 
tion within 1 R e (see ?), even though NGC 4486 (M 87) is in fact 
flattened at very large radii (see ?, and references therein). This 
implies that only a few galaxies near the Xr c — 0.1 threshold 
(NGC 5982, NGC 4550, NGC 4278) could cross that threshold if 
we were to have a complete coverage up to 1 Re ■ 

This is also illustrated in Fig. [3] which shows that histograms 
of Xr c for radii of 1 R e and R e /2 (or restricted to the equiva- 
lent effective aperture of the SAURON field of view, whichever is 
smaller) have the same fraction of slow and fast rotators. The over- 
all distribution of An values is similar for R e and R e /2, although 
Xr tends to increase significantly from R e /2 to R e for fast rotators 
(see also the solid straight lines in Fig. [5j, an obvious implication 
of the observed rising An profiles in Fig. [2] There are 36 fast rota- 
tors and 12 slow rotators (75% and 25% of the total sample), their 
median Xr c values being ~ 0.44 and 0.05 respectively. Within 
the class of slow rotators, three galaxies have Xr c significantly be- 
low 0.03 (their mean stellar velocity maps being consistent with 
zero rotation everywhere, as mentioned above). These are among 
the brightest galaxies of our sample, namely NGC 4486, NGC 4374 
and NGC 5846. 

3.3 Misalignments and twists 

The fact that slow and fast rotators exhibit distinct kinematics can 
be demonstrated by considering the global alignment (or misalign- 
ment) between photometry and kinemetry. Fig.|4]illustrates this by 
showing the kinematic misalignment ^ = \PA p hot — PAti n \ for 
all 48 E and SO galaxies in the SAURON sample. Note that the 
photometric PA is derived using the large-scale MDM data, but 
the kinematic PA via a global measurement on the SAURON ve- 
locity maps as described in Appendix C of ?. All fast rotators, ex- 
cept one, have misalignments below 10°. The only exception is 
NGC 474 which is an interacting galaxy with well-know irregular 
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Figure 4. Xr s versus the kinematic misalignment between the global 
photometric major-axis and the kinematic axis within the SAURON field. 
Slow rotators are represented by red circles, fast rotators by special blue 
symbols (horizontal ellipse plus a vertical line). The vertical dashed line 
corresponds to ^ = 5°. Nearly all fast rotators have small * values 
(< 10°), the only exception being NGC474, the photometry of which is 
perturbed by the presence of irregular shells. This contrast with slow rota- 
tors which show significantly non-zero values. 

shells (?). In fact, the few galaxies which have 5° < * < 10° 
(NGC 3377, NGC 3384, NGC 4382, NGC 4477, NGC 7332) are al- 
most certainly barred. Even the two relatively face-on galaxies 
NGC 4262 and NGC 4477, the photometry of which shows the sig- 
nature of a strong bar within the SAURON field of view, have their 
outer photometric PA well aligned with the measured kinematic 
PA. In contrast, more than half of all slow rotators have ^ > 10°, 
and none of these exhibit any hint of a bar. This difference in the 
misalignment values of slow and fast rotators cannot be entirely due 
to the effect of inclination, not only because of the argument men- 
tioned in Sect. 14.21 but also because even the roundest fast rotators 
do not exhibit large misalignment values (see Sect. 5. 1 in Paper X). 

Another remarkable feature comes from the observed veloc- 
ity twists in the SAURON kinematic maps: only 6 galaxies out of 
48 exhibit strong velocity twists larger than 30° outside the inner 
3", namely NGC 3414, 3608, 4550, 4552, 5198, and 5982, with 
3 out of these 6 having large-scale counter-rotating stellar compo- 
nents (NGC 3414, 3608, and 4550). All these galaxies are in fact 
slow rotators. This implies that only fast rotators have a relatively 
well defined (apparent) kinematic major-axis, which in addition is 
roughly aligned with the photometric major-axis. 

There is a close link between Xr and the (apparent) specific 
angular momentum (see AppendixfA}- Xr is therefore a continuous 
parameter which provides a quantitative assessment of the apparent 
mean stellar rotation. As emphasised in this Section, slow and fast 
rotators exhibit quantitatively but also qualitatively different stellar 
kinematics. This strongly suggests that slow rotators as a class can- 
not simply be scaled-down versions of galaxies with Xr c > 0.1. 

3.4 Inclination effects 

Obviously, Xr is derived from projected quantities and therefore 
also significantly depends on the viewing angles. Can slow rotators 



be face-on versions of fast rotators? Assuming Xr roughly follows 
the behaviour of V/a with the inclination angle i (Appendix IB!, a 
galaxy with a measured Xr f = 0.05, typical of observed slow rota- 
tors, would require to be at a nearly face-on inclination of i ~ 20° 
to reach an intrinsic (edge-on) value of Xr c ~ 0.15, the smallest 
value for all observed fast rotators in our sample. The probability 
that slow rotators are truly fast rotators seen face-on is therefore 
small, and cannot explain the 25% population of slow rotators ob- 
served in the Xr f histogram (Fig. [3}- 

As emphasised above, the measured Xr c values are most 
probably lower limits because we observe galaxies away from 
edge-on. An illustrative example is provided by the case of 
NGC 524. According to Paper IV, NGC 524 is viewed with an in- 
clination of 19°: still, NGC 524 exhibits a regular velocity pattern, 
and is a fast rotator with An e ~ 0.28. There is also some evidence 
that the three fast rotators with the lowest Xr c values in our sample, 
namely NGC 3379, NGC 4278 and NGC 4382, are significantly in- 
clined galaxies (?; Papers III, IV and ?, hereafter Paper V; ?). This 
illustrates the fact that the apparent distribution of Xr c within the 
fast rotator class may not fully reflect the intrinsic distribution with 
galaxies viewed edge-on. It is, however, impossible to estimate the 
"edge-on" Xr s distribution for fast rotators as this would require 
ad minima an accurate measurement of their inclination. Although 
the fast rotator class may probably include truly mild as well as 
very fast rotators (as expected from the continuous nature of A_r c ), 
all galaxies with Xr c > 0.1 exhibit a global and regular rotation 
pattern, with a clear sense of rotation and a significant amount of 
specific stellar angular momentum. This, to us, justifies the use of 
the label "fast rotators" for these galaxies. Further discussions of 
the true nature of fast rotators as a class has to wait for a larger and 
complete sample. 



4 OTHER CLASSIFICATIONS 

In this Section, we examine in more detail other galaxy properties, 
including photometry and kinemetry, to determine if these could 
help in defining kinematic classes, and their potential link with Xr. 

4.1 Hubble Classification 

There is a clear overlap of Es and SOs in the fast rotator class, with 
most of the corresponding galaxies having ellipticities higher than 
0.2. This illustrates the fact that many Es have kinematic charac- 
teristics similar to SOs (??). Some galaxies classified as ellipticals 
also show photometric signatures of embedded discs, quantified as 
a positive a&/a, and were classified as "discy" ellipticals or E(d) 
by KB96. It seems difficult, however, to distinguish between a so- 
called E(d) and an SO galaxy from the photometric properties or 
the anisotropy diagram alone. KB96 argued that E(d)s are objects 
intermediate between SOs and boxy ellipticals or E(b)s (ellipticals 
with boxy isophotes, i.e. negative a^/a). In principle, the Ed-SO- 
Sa sequence is one of decreasing bulge-to-total light ratio (B/T), 
with SOs having B/T ~ 0.6 in the B band (?). B/T is, however, a 
fairly difficult quantity to measure in early-type galaxies, as it de- 
pends on the adopted model for the surface brightness distribution 
of the disc and bulge components. As emphasised by ?, the kine- 
matic information is critical in assessing the rotational support of 
both components. 

Consider two galaxies of our sample, NGC 3377 and 
NGC 2549, which both have similar Xr p , the former being clas- 
sified as a discy elliptical, the latter as a lenticular. Both have very 
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Table 1. Characteristics of the E and SO galaxies in the representative SAURON sample. All galaxies with A^ e > 0.1 are classified as fast rotators. 
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Notes: (1) Galaxy identifier (NGC number). (2) Hubble type (NED). (3) Numerical morphological type (LEDA ?). (4) Galaxy distance modulus from ??, 
(corrected by subtracting 0.06 mag, see ?), or from the LEDA database assuming a Hubble flow with H = 75 kms" 1 Mpc" 1 . (5) Absolute B magnitude 
(Paper II). (6) Effective radius, in arcsec. (7) Velocity dispersion derived using the luminosity-weighted spectrum within R e or within the SAURON field, 
whichever is smaller. (8) and (9) Global ellipticity within R e /2 and R e (from SAURON), or within the SAURON field, whichever is smaller. (10) Mean 
ellipticity within 1 R e , derived from the GALPHOT radial profile. (11) Mean isophote shape parameter a±/a (in %) within 1 R e . (12) and (13) V/cr (see 
Paper X) and \r within 1 R e . (14) Kinemetry group (see Sect. 12.41 . (15) Rotator class: F=fast, S=slow. 



similar total luminosity, V — I colour (?) and gas content (Paper V). 
Their SAURON stellar kinematic maps are also quite similar (Pa- 
per III), with large-scale disc-like rotation, a centrally peaked stellar 
velocity dispersion and a significant /13 term, anti-correlated with 
the mean velocity. Finally, there is evidence that NGC 3377 con- 



tains a bar, with its velocity map exhibiting a stellar kinematical 
misalignment and a spiral-like ionised gas distribution (Paper V). 
This therefore strongly argues for NGC 3377 to be a misclassified 
barred SO (SAB0). Many authors (e.g. ???, and KB96) indeed sug- 
gested that many Es are in fact misclassified SOs, mostly on the 
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Figure 5. \r c versus the global ellipticity e for the 48 E and SO galaxies 
of the SAURON sample. The solid lines indicate where the points would 
move for a smaller equivalent effective aperture of R e /2. Slow rotators 
(which have A/j e ^ 0.1) are red circles, and fast rotators (A^ c > 0.1) are 
represented by special blue symbols (horizontal ellipse plus a vertical line). 
The black dashed line corresponds to the curve expected for isotropic oblate 
rotators viewed edge-on (see Appendix 151 and ?). The galaxy NGC4550 is 
labelled. 



basis of photometry alone. Considering the kinematic properties of 
the observed SAURON early-type galaxies, we suspect that most 
and possibly all of the 13 Es which are fast rotators are in fact mis- 
classified SO galaxies, a conclusion also supported by the results of 
Paper X via the use of state-of-the-art dynamical models. 
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Figure 6. Mean ellipticity e e (top) and Xn r (bottom) versus mean (0.4, /a) e 
(in %) within 1 R e . Symbols for slow and fast rotators are as in Fig. [5] The 
vertical dashed lines correspond to (04/0) e = ±1% 



4.2 Isophote shapes 

Figure [5] shows Xr c versus the global ellipticity e for the 48 
SAURON E and SOs, measured using an aperture of 1 R e (or in- 
cluding the full SAURON field of view for galaxies with large R e , 
see Sect. l3.2b . All galaxies in Fig.[5]lie close to or below the curve 
expected for isotropic oblate rotators viewed edge-on (see Ap- 
pendix|B]and ?). Fast rotators, which have velocity maps with sig- 
nificant large-scale rotation, have ellipticities ranging up to about 
0.6. Apart from the atypical case of NGC 4550, slow rotators show 
little or spatially confined rotation and all have ellipticities e < 0.3. 
NGC 4550 is in fact a nearly edge-on galaxy with two co-spatial 
counter-rotating stellar discs, each contributing for about 50% of 
the total luminosity, and should therefore be regarded as an atypi- 
cal case where a high ellipticity is accompanied by a relatively low 
mean stellar velocity in the equatorial plane (??, and see Paper X 
for a detailed discussion). 

KB96 mention that 0,4, /a can be used as a reasonably reli- 
able way to measure velocity anisotropy. They first observe, us- 
ing (V/cr)* (V/cr normalised to the value expected for an isotropic 
edge-on oblate system), that discy galaxies (0,4, /a > 0) seem con- 
sistent with near isotropy, whereas boxy galaxies (04/ a < 0) 
spread over a large range of (V/cr)* values. However, as empha- 
sised by ? and in Paper X, (V/cr)* is not a good indicator of 
anisotropy: the relation between (V/cr)* and the anisotropy of a 
galaxy strongly depends on its flattening, as well as on its incli- 
nation (?). More flattened galaxies will thus tend to lie closer to 



the (V/cr)* — 1 curve, which explains in part why galaxies with 
discy and boxy isophotes seem to extend over different ranges 
of (V/cr)*, the former being then incorrectly interpreted as near- 
isotropic systems. This suggests we should avoid using (V/cr)* at 
all, and leads us to examine more directly the orbital anisotropy 
of galaxies (see ?). This is achieved in Paper X, in which a trend 
between the anisotropy parameter in the meridional plane and the 
average intrinsic ellipticity is revealed (see their Fig. 7), intrinsi- 
cally more flattened galaxies tending to be more anisotropic. This 
contradicts the view that discy galaxies are nearly isotropic, and 
we should therefore reexamine the relation between the isophotal 
shape as measured by 04/ a and the kinematic status of early-type 
galaxies. 

KB96 presented a fairly significant correlation between the 
mean ellipticity e and 04/ a for their sample of early-type galaxies 
(their Fig. 3). They mentioned that galaxies in a e versus (14/0 dia- 
gram exhibit a V-shaped distribution. The authors therefore pointed 
out that this remarkable correlation implies that early-type galaxies 
with nearly elliptical isophotes are also nearly round, and galaxies 
with the most extreme a4/a values are also the most flattened (see 
also ?). We confirm this result using our photometric data on the 
SAURON sample as shown in Fig. [6] (upper panel). As mentioned 
above, all slow rotators with the exception of NGC 4550 have el- 
lipticities e e lower than 0.3, and our expectation that these galaxies 
should have isophotes close to pure ellipses is verified: 10 out of 1 1 
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Figure 7. A^ e versus absolute magnitude Mb (left panel) and virial mass M v i r (right panel) for the 48 E and SO of the SAURON sample. In the left panel, 
symbols correspond to the cusp slope classification (????) with power-laws as triangles, cores as squares, NGC 821 which is an "intermediate" object as a 
circle, and crosses indicating galaxies for which there is no published classification. In the right panel, symbols for slow and fast rotators are as in Fig. [5] and 
NGC numbers for a few galaxies are indicated. 



have |(a 4 /a)e| < 1% (with 9 out of 1 1 having J(a 4 /a) e | < 0.5%), 
the exception being NGC 3414, for which the relatively large posi- 
tive (0,4/(1),, value comes from the remarkable polar-ring structure 
(?)■ 

In fact, most fast rotators have a maximum absolute | (a 4 /a) | 
value within 1 R e larger than 2%, whereas most slow rotators (with 
the exceptions again of NGC 3414 and NGC 4550) have maximum 
absolute |(a4/a)| values less than 2%. This means that most fast 
rotators in the SAURON sample of E and SO galaxies exhibit sig- 
nificantly non-elliptical isophotes, but most slow rotators do not. 
There is a tendency for the fast rotators to have positive 04 /a (discy 
isophotes), and most slow rotators (7 out of 12) have negative 0,4 /a 
(boxy isophotes) We do not detect any global and significant cor- 
relation between the boxiness and the angular momentum per unit 
mass as quantified by Xr c in the galaxies of our sample, so that 
I (0,4/ a) I is clearly not a good proxy for rotation or angular momen- 
tum. This seems partly in contradiction with the claim by KB96 that 
rotation is dynamically less important in boxy than in discy early- 
type galaxies. The latter result probably originated from the com- 
bination of two observed facts: firstly boxy galaxies are on aver- 
age less flattened (extremely flattened galaxies are very discy), and 
secondly the use of (V/a)* as a proxy for the amount of organised 
rotation in the stellar component, which tends to underestimate the 
rotational support for less flattened galaxies. 



4.3 Luminosity and mass 

In the left panel of Fig.|7] we show the distribution of A.r c as a func- 
tion of absolute magnitude M B for the 48 SAURON E/S0 galax- 
ies. The three slowest rotators (NGC 4486, NGC 4374, NGC 5846) 
are among the brightest galaxies in our sample with Mb < 
—21 mag. Other slow rotators tend to be bright but are spread over 



a wide range of absolute magnitude going from the relatively faint 
NGC 4458, to brighter objects such as NGC 5813. The bright and 
faint end of slow rotators can be distinguished by the shapes of 
their isophotes: slow rotators brighter than Mb < —20 mag all 
exhibit mildly boxy isophotes (negative 04 /a with amplitude less 
than 1%) while the four discy slow rotators are all fainter than 
Mb > —20 mag, following the known correlation between the 
isophote shapes and the total luminosity of early-type galaxies (?). 
Most fast rotators are fainter than Mb > —20.5 mag. 

Going from total luminosity to mass, we have estimated the 
latter by approximating it with the virial mass M„j r derived from 
the best-fitting relation obtained in Paper IV, namely M V i r ~ 
5.0 R e a1/G, where a e is the luminosity weighted second veloc- 
ity moment within 1 R e (see Paper IV for details). The calculation 
of M V i r from observables depends on the distance of the object, 
which we obtained from different sources for the galaxies in our 
sample (in order of priority, from ??, and from the LEDA database 
assuming a Hubble flow with H = 75 km s _1 Mpc -1 ). A trend of 
Xr c tending to be lower for more massive galaxies clearly emerges 
if we now use this estimate of the virial mass M V i r (right panel of 
Fig. 0, as expected from the one observed with absolute magni- 
tude Mb (left panel of Fig.[7j>. The three slowest rotators are in the 
high range of M V i r with values above lO 11 ' 5 Mq. There is a clear 
overlap in mass between fast and slow rotators for M„; r between 
10 11 and 10 11 ' 5 M . However, all slow rotators, besides NGC 4458 
and NGC 4550, have M vir > 10 11 M , whereas most fast rotators 
have M V i r < 10 Mq, lower masses being reached as the value of 
Xr c increases. It is worth pointing out that the absolute magnitude 
in the K band, Mk , very nicely correlates with M„, r (significantly 
better than with Mb), so that the trend observed between Xr and 
M v i r is also present if we were to examine the relation between Xr 
and Mk ■ 

Out of the 48 SAURON galaxies, 33 have published cusp 
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Figure 8. \r versus mean ellipticity e, with indications of the velocity 
structures identified in Sect. 12.41 Galaxies with MCs are shown as dia- 
monds, SCs as circles, and CLVs as right triangles. Colours for slow and 
fast rotators are as in Fig. \5\ KDC stand for kinematically decoupled com- 
ponents, KM and KT for kinematic misalignment and twist, respectively, 
LV for Low Velocity component. 

slope classification, distinguishing "core" and "power-law" galax- 
ies (????). As illustrated in the left panel of Fig.[7] we find that most 
slow rotators are core galaxies, and most fast rotators are power- 
law galaxies (note that NGC 821 has an A_R e ~ 0.25 and is spec- 
ified as an "intermediate" object between a core and power-law in 
?). There are yet no core galaxies with A_r c > 0.3, although the 
inclined galaxy NGC 524 (see Paper IV) would have a very high 
A_R e value if seen edge-on. These results are expected as there is 
a known trend between the central luminosity gradient and the to- 
tal luminosity of early-type galaxies, bright members tending to be 
core galaxies, and lower luminosity ones to have power-law profiles 
(?). All galaxies with X Re > 0.3 have indeed Mb > -20.7 and 
conversely none have M V i r > 10 Mq. This is, however, not 
a one-to-one correspondence since we find both power-law in slow 
rotators (NGC 3414 and NGC 5813), and core-like fast rotator (e.g., 
NGC 524). There is also a domain in mass and luminosity where 
both cusp and cores are found, namely for M v i r between 10 and 
10 11 ' 5 M Q , or correspondingly for M B from -19.5 to -20.5. An- 
other interesting result is found when examining the larger-scale 
luminosity profiles of galaxies in our sample via the representa- 
tion by a Sersic law. Besides the atypical case of NGC 4550, all 
slow rotators have Sersic index n > 4: again, this is expected since 
galaxies with larger Sersic shape index tend to be brighter (??). Fi- 
nally, galaxies with the lowest Sersic n values are also among the 
fastest rotators. A more detailed account regarding these issues will 
be provided in Falcon-Barroso et al. (in preparation). 

4.4 Kinemetry groups 

We now turn to the kinemetric profiles (see Sect. I2.4t of the 48 
SAURON galaxies, which allow us to determine the number of ob- 
served kinematic components and their individual characteristics. 
We thus make use of the average photometric and kinematic posi- 
tion angles and axis ratio (PA p hot and PAfej n , q and qum), and the 
average and maximum of the velocity amplitude fei. 

The kinemetric groups of the observed 48 SAURON velocity 



maps, defined in Sect. 12.41 are provided in Table[T] Out of 48 galax- 
ies in this sample, 33 (69%) exhibit multiple components (MC), in- 
cluding 10 kinematically decoupled components (KDCs; 21%) and 
4 objects with central low-level velocity (CLVs; namely NGC 3032, 
NGC 4150, NGC 4382, and NGC 7457). Two more galaxies with 
KDCs are in fact at the limit of being CLVs (NGC 4621 and 
NGC 7332) with a maximum velocity for their inner component 
around 20 km s~ at the SAURON resolution. The latter two galax- 
ies as well as the CLVs have in fact all been shown to harbour small 
counter-rotating stellar systems (???, hereafter Paper VIII): it is the 
lower spatial resolution of the SAURON data which produces the 
low central velocity gradient. If we therefore count CLVs as galax- 
ies with KDCs, this would result in a total of 14 KDC (29%) early- 
type objects. Among the KDCs, 5 galaxies (NGC 3414, NGC 3608, 
NGC 4458, NGC5813, NGC5831) have outer LV components. 
Kinematic misalignments are observed in 12 galaxies (25%), and 
6 galaxies have individual kinematic components which have kine- 
matical twists (KTs, 2 of them also having kinematical misalign- 
ments - KM). Considering the difficulty of detecting these struc- 
tures at certain viewing angles or low spatial resolution, the number 
of such detected velocity structures is a lower limit. 

In Fig. [8] we show the (A_R e ,e) diagram (as in Fig. [5}, now 
with the characteristic velocity structures detected via kinemetry. 
The three slowest rotators are obviously tagged as low-velocity 
(LV) systems. As mentioned above, these three galaxies are gi- 
ant (bright and massive) roundish ellipticals with Mb < —21. 
Apart from these three and the atypical case of NGC 4550, all 
other slow rotators harbour KDCs, the central kinematic compo- 
nent having a typical size of 1 kpc or larger (see Paper VIII). 
This contrasts with KDCs and CLVs in fast rotators (NGC 3032, 
NGC 4150, NGC 7332, NGC 7457), where the size of the central 
(counter-rotating) stellar component is only a few arcseconds in ra- 
dius, corresponding to significantly less than 500 pc (Paper VIII). 

A number of fast rotators show kinematic misalignments 
(KM) and twists (KT), most of these being relatively small in 
contrast with the ones found in slow rotators (see Sect. 13. 3\ . The 
four fastest rotators with a KT, KM, or KDC, namely NGC 1023, 
NGC 3377, NGC 7332 and NGC 7457, are all barred galaxies. Most 
of the galaxies with multi-components but no specific velocity 
structures are among the fastest rotators. This result may be partly 
understood if these objects tend to be close to edge-on, which ren- 
ders the detection of such velocity structures harder, but the detec- 
tion of multiple components easier. 



5 DISCUSSION 

The fast rotators are mostly discy galaxies exhibiting multiple com- 
ponents in their stellar velocity fields. Contrarily to slow rotators, 
the main stellar kinematic axis in fast rotators is relatively well 
aligned with the photometric major-axis, except for one galaxy 
(NGC 474) known to harbour irregular shells. This result together 
with the fact that the \r profiles are qualitatively different for slow 
and fast rotators clearly show that slow rotators are not velocity 
scaled-down versions of fast rotators. The small number of galax- 
ies in our sample, as well as our biased representation of the galaxy 
luminosity function, reminds us that a larger and complete sample 
is required to reveal the true A_R e distribution in early-type galaxies. 

We therefore first briefly discuss here the results obtained so 
far, this time including 18 additional early-type galaxies ("spe- 
cials") observed with SAURON within the course of other specific 
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Figure 9. Left panel: A^j e versus the ellipticity e including the additional 1 8 
E and SO galaxies observed with SAURON (labelled with NGC numbers). 
Symbols and the dashed line are as in Fig. [5] Right panel: histogram of A^ e 
values including the 18 "specials". 



projects. For these additional targets, we had to retrieve the main 
photometric parameters (e.g., R e , e, qa/cl) from the literature (e.g. 
??), the kinematic measurements (cr e , A_R e ) being derived from the 
available SAURON data as for the 48 E/SO galaxies of the main sam- 
ple (only 12 galaxies out of the 18 specials have available 04 mea- 
surements). Adding these 18 galaxies does not change the overall 
distribution of fast and slow rotators in a A_r e versus e diagram, as 
shown in Fig. [9] The fraction of slow rotators (17 out of 66) is still 
~ 25%. We also confirm that most slow rotators have relatively 
small ellipticities (e < 0.3), in contrast with fast rotators which 
span a wide range of flattening. The most remarkable fact lies in the 
confirmation that within these 18 extra targets, there are 5 slow ro- 
tators (NGC 4168, NGC 4406, NGC 4472, NGC 4261, NGC 4365), 
and they all contain large (kpc) scale KDCs. Similarly, in Fig. [10] 
we confirm the lack of a global significant correlation between \n s 
and a^/a as shown in Fig. [6] Again, about two thirds of the slow 
rotators are boxy (5 out of 17 are discy, and 1 has 0,4 /a — ), and 
very discy galaxies tend to be fast rotators with high A_R e values. 
However, 04 /a is clearly not a good proxy for angular momentum. 
We then show in Fig.[TT]the relation between Xjt e and M V i r with an 
additional 17 "specials" (excluding the compact galaxy NGC 221 
[M 32] which would stretch the plot unnecessarily). Most slow ro- 
tators are massive galaxies with M V i r > 10 11 Mq, with still only 
two exceptions (NGC 4550 and NGC 4458), reinforcing the results 
from Fig. QT] The trend for more massive galaxies to have lower 
A_b c values is still observed, with an overlap of fast and slow rota- 
tors in the range 10 11 - 10 12 M Q . The fast rotator NGC 2320 seems 
to have a rather large A.r c for its virial mass (being in fact the most 
massive galaxy out of the 66). NGC 2320 is the most distant galaxy 
of this set (with D ~ 83 Mpc), and has a rather unusual molecu- 
lar gas content for an early-type galaxy: it exhibits an asymmetric 
molecular gas disc with a mass of about 4 x 10 9 Mq, interpreted 
as a sign of recent accretion or dynamical perturbation (?). 

In Fig.Qj] we also show the "cusp / core" classification as in 
Sect. 14.31 Most fast rotators are still galaxies with cusps, and most 
slow rotators are core galaxies. Conversely, massive galaxies tend 
to have cores, and smaller ones tend to exhibit cuspy central lumi- 
nosity profiles: in fact all galaxies with M„j r > 10 11 ' 5 Mq have 
cores. There is a transition region, in the range 10 11 — 10 11 ' 5 Mq 
where we find both cusp and core galaxies. More interestingly, 
all core galaxies still have Aj? e < 0.3. However, a galaxy like 
NGC 524 is almost certainly very inclined (see Paper IV), its cor- 
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Figure 10. Xjt e versus 04 for the 48 E/S0 of the SAURON sample and an 
additional 12 galaxies for which we have available 04 values. Symbols for 
slow and fast rotators are as in Fig. [6] 
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Figure 11. \ Re versus the virial mass M vir for the 48 E/SO of the SAURON 
sample plus an additional 17 "specials" (we do not include the fast rotator 
NGC 221 - M 32 - in this plot considering its very low mass). The top panel 
shows histograms of M v i r (in log, with steps of 0.5) for both slow (red) and 
fast rotators (blue). Colours in both panels for slow and fast rotators are as 
in Fig. [6] NGC numbers are indicated for all specials. Symbols for galaxies 
with cores and cusps are as in the left panel of Fig. [7] 



responding "edge-on" A_r s value being then much larger than 0.75. 
It would thus be interesting to understand if an edge-on version of 
NGC 524 would still have a core-like central luminosity profile. We 
also observe a few galaxies with cusps and rather low \r c values 
(NGC 5831) which cannot be reconciled with inclined fast rota- 
tors. We may therefore only discriminate cusp and core galaxies 
combining both their mass and their Xr c values: massive galaxies 
(M vtr > 10 115 Mq) with A flc < 0.3 are indeed all core galaxies, 
while smaller galaxies (M„i r < 10 11 Mq) with A_b c > 0.3 seem 
to all be galaxies with cusps. The fact that we observe very few 
massive (M„; r > 10 11 Mq) galaxies with high Au c values implies 
that either our sample is biased against them (e.g. we only observed 
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close to face-on massive galaxies), or these galaxies are rare. Since 
it is difficult to understand how such a bias could affect our sample, 
we need to conclude that indeed massive galaxies tend to have little 
or no baryonic angular momentum within one effective radius, and 
that the trend observed in Fig.QT]is real (see e.g., the upper panel 
with the histogram of M„ir)- 

Interpreting the measured kinematic parameter Xr as a proxy 
for the amount of stellar angular momentum per unit mass in the 
central region of early-type galaxies (see Appendix A), we can then 
discuss the origin of this angular momentum. The standard scenario 
for the formation of galaxy structures includes hierarchical cluster- 
ing of cold dark matter halos within which gas is cooling (???). The 
angular momentum of the dark matter halos is thought to originate 
in cosmological torques and major mergers (e.g. ?), this hypothe- 
sis providing a reasonable frame for the formation of disc galaxies. 
If major mergers produce a significant increase in the specific an- 
gular momentum of the dark matter halos at large radii (?), minor 
mergers seem to just preserve or only slightly increase it with time 
(?). But little is known on the expected distribution of the angular 
momentum of the baryons (??), and even less if we focus on the 
central regions of galaxies (within a few effective radii). In fact, 
discs formed in numerical simulations are generally an order of 
magnitude too small (but see ?, for a possible solution). 

Numerical simulations, whether they include a dissipative 
component or not, have helped us to understand how mergers influ- 
ence the rotational support of the baryonic component in galaxies 
(see e.g. ?????). Stellar discs are cold and fragile systems, thus eas- 
ily destroyed during major mergers which often lead to elliptical- 
like remnants even in the presence of a moderate amount of gas 
(?). Intermediate to minor mergers preserve part of the disc better, 
but in most cases a merger leads to a redistribution of the angular 
momentum of the central stellar component outwards (?). 

In this context, fast rotators have either preserved or regained 
their specific angular momentum in the central part. Since both 
very gas-rich major and minor mergers seem to produce fast (disc- 
dominated) rotators (??), this requires either the absence of a major 
dry merger which would expel most of the angular momentum out- 
wards, or the rebuilding of a disc-like system via gas accretion. At 
high redshift (z > 2), gas was abundant, and very gas-rich mergers 
should have therefore been common. ? in fact claim that progeni- 
tors of early-type galaxies must be gas-rich (gas fraction > 30%) to 
produce the tilt of the fundamental plane (FP). Dry major mergers 
are therefore expected to occur preferentially at lower z, with fast 
rotators not having suffered from such events. In the picture of the 
hierarchical formation of structures, minor mergers are more com- 
mon than major mergers and we can thus expect galaxies (slow and 
fast rotators) to have suffered from more than one of these events 
up to z = 0. 

The resulting specific angular momentum of fast rotators, as 
quantified by Xr, is therefore expected to result mostly from a com- 
petition between (i) gas-rich minor mergers or other inflow of ex- 
ternal gas that causes a gradual increase of A_r (in the inner parts), 
and (ii) dissipationless dry minor mergers triggering disc instabil- 
ity and heating, and resulting in the transformation of disc material 
into a more spheroidal component, lowering Xr (in the inner parts). 
The scenario described here follows the idea previously sketched 
by many authors (e.g. ????) that dissipation is important in the for- 
mation process of fast rotating early-type galaxies. The 'heating' of 
the disc via star formation leads to an increase in the vertical dis- 
persion, more isotropic and rounder galaxies, and therefore mov- 
ing the galaxies along the anisotropy-ellipticity trend (Paper X) 
towards the slow rotators. A sudden removal of the gas (e.g., due to 
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Figure 12. Average luminosity-weighted age of the stellar component ver- 
sus the specific mass of ionised gas within the SAURON field of view 
Mgas/Mui,. (ionised gas mass normalised by the virial mass, see Sect. 14.31 . 
Symbols for fast and slow rotators are as in Fig. [3] Numbers for the 
ionised gas specific mass are derived from values in Papers IV and Pa- 
per V, and average luminosity-weighted ages are from Kuntschner et al. 
(in preparation). The vertical and horizontal dashed lines shows the limits 
of > 10 -6 MgasfM vir and 6.5 Gyr (z ~ 1), respectively. 



AGN feedback or ram pressure stripping) might have 'quenched' 
this process, and quickly moved the galaxy from the 'blue cloud' 
to the 'red sequence' (?). Some galaxies may still have recently 
accreted some gas (e.g. via interaction with a companion, see ?), 
and sometimes show the presence of a younger stellar population 
(such as NGC3032, NGC3489, NGC4150, see ?, Paper VI, and 
see also Paper VIII). In fact all SAURON E/SO galaxies which have 
luminosity weighted ages of their stellar component lower than 
6.5 Gyr (z ~ 1) are fast rotators (Fig. [12] see also Kuntschner 
et al. in preparation). Among the 11 galaxies in our sample with 
a high ionised gas content (M gas /M V i r > 10~ 6 ), only one is a 
slow rotator, namely the atypical disc galaxy NGC 4550. The slow 
rotator with the second highest gas content is the polar-ring galaxy 
NGC 3414 (M gas /M vir ~ 0.6 x 10~ 6 ), NGC 4550 and NGC 3414 
being in fact the two galaxies having deviant (a4/a) e values with 
respect to all other slow rotators (see Sect. |4.2t . 

Following the same line of argument, slow rotators need to 
have expelled most of their angular momentum within one effec- 
tive radius outwards. This would require a significant fraction of 
the mass to be accreted during mergers with relatively gas-poor 
content. As emphasised above, all slow rotators with A_r c > 0.03 
(thus excluding the three slowest rotators which are consistent with 
zero rotation) have stellar kpc-size KDCs. Opposite to the smaller- 
scale KDCs in fast rotators (Sect. pOt , the latter have a similar old 
stellar population as the rest of the galaxy (Paper VIII), and thus 
were formed long ago when gas was more abundant than today. 
Indeed, numerical simulations of the formation of such structures 
seem to require the presence of some gas during the merger pro- 
cess followed by subsequent star formation, either in equal-mass 
mergers (?) or in multiple intermediate or minor mergers (Bour- 
naud, PhD Thesis, Paris). Fig. [12] reveals that most slow rotators 
have a relatively modest fraction of ionised gas and are in fact 
older than ~ 6.5 Gyr (z ~ 1). The only two exceptions are, again, 
NGC 4550 which is very probably the result of the rare encounter 
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of two disc galaxies with nearly exactly opposite spins (see Pa- 
per X), and NGC 3414 which shows the "rudiments of a disc" (?), 
and photometric structures reminiscent of polar-ring galaxies (?) 
presumably formed via a tidal gas accretion event (?). These re- 
sults are consistent with the suggestion that slow rotators did not 
suffer from significant (in mass) recent dry mergers which may de- 
stroy or diffuse the large-scale KDCs. Minor mergers may still have 
recently occurred, and as above result in a competition of (slightly) 
lowering and increasing Xr. 

The three slowest rotators (NGC 4374, NGC 4486, and 
NGC 5846) may then correspond to the extreme end point of early- 
type galaxy evolution where dry mergers had a major role in the 
evolution, allowing Xr to reach values consistent with zero in the 
central region. These galaxies are found in rather dense galactic en- 
vironment and are expected to have experienced a significant num- 
ber of mergers with galaxies which were themselves gas-depleted 
(due to ram-pressure stripping, due to efficient AGN feedback, or 
due to earlier mergers), preferentially on radial orbits along the cos- 
mological filaments (?). As for the other slow-rotators, the initial 
mergers will have been gas-rich to place them on the fundamen- 
tal plane, but later-on dry mergers are expected to be frequent as 
galaxies fall into the cluster potential well. The insignificant role 
of star formation in the late history of these galaxies is reflected 
in both their low specific gas content and old luminosity weighted 
ages (Fig.ll2t. Overall, the trend we observe for early-type galax- 
ies to have lower A r with increasing mass (Fig. O is consistent 
with the gradually more important role of (gas-poor) mergers when 
going from fast rotators, to slow rotators. 

Is this also consistent with the observed trend between Xr^, 
M V i r and the cusp slope described above? Cores are thought to 
be the result of the scouring via a binary supermassive black hole 
(??, and references therein), the potential consequence of a merger. 
After a core is formed, a subsequent central accretion of gaseous 
material (from external or internal sources) followed by an episode 
of star formation could rejuvenate a cuspy luminosity profile, most 
probably in the form of an inner disc. If this scenario is correct, 
there should be a relation between the presence of a core and Xr c . 
The observation that massive slow rotators, thought to have signif- 
icantly suffered from dissipationless mergers, all exhibit cores, and 
that small galaxies with high Aij c all contain cusps does indeed fit 
into this scenario. However, the fact that minor gas-free mergers 
may only remove part of the baryonic angular momentum, and that 
gas can be accreted from e.g., external sources or stellar mass loss, 
implies that there should not be a one-to-one relation, as indeed ob- 
served. Intriguingly, the two slow rotators with cusps (NGC 3414, 
NGC 5813) both show the presence of a decoupled stellar disc-like 
component in their central few hundreds of parsecs (as witnessed 
by the amplitude of the associated /t3 parameter, see Paper III). 
These may illustrate again the competition between gas-poor merg- 
ers and gas accreting processes in the making of galaxy structures. 



and quantitatively different stellar kinematics (see also Paper X). 
Galaxies with A.r ^0.1 form a distinct class characterised by lit- 
tle or no large-scale rotation, by the presence of kinematically de- 
coupled cores, and by significant kinematic misalignments and/or 
velocity twists. This contrasts with galaxies for which A_r > 0.1, 
which display global rotation along a well defined apparent kine- 
matic major-axis that is nearly aligned with the photometric major- 
axis. We have also shown that slow and fast rotators exhibit qualita- 
tively different An radial profiles and that slow-rotators are unlikely 
to be face-on versions of fast-rotators. 

Beside the obvious difference in their angular momentum con- 
tent, the main properties of the fast and slow rotators are as follows: 

(i) 75% of the galaxies in our sample are fast rotators, 25% are 
slow rotators. Most of the slow rotators are classified as Es, while 
fast rotators include a mix of Es and SOs. 

(ii) Within the slow rotators, three early-type galaxies have ve- 
locity fields consistent with near zero rotation. These three objects 
(NGC4374=M84, NGC4486=M87 and NGC 5846) are bright 
and massive, nearly round galaxies. 

(iii) Fast rotators have ellipticities up to e ~ 0.6, in contrast to 
slow rotators which are relatively round systems with e < 0.3. 

(iv) Fast rotators tend to be relatively low-luminosity objects 
with Mb > —20.7. Slow rotators cover almost the full range of 
absolute magnitude of our sample, from faint discy galaxies such 
as NGC 4458 to bright slightly boxy galaxies such as NGC 5813, 
although they on average tend to be brighter than fast rotators. All 
slow rotators have Sersic index n larger than 4, except the atypical 
galaxy NGC 4550. 

(v) More than 70% of the galaxies in our sample exhibit multi- 
component kinematic systems. As mentioned above, fast rotators 
have nearly aligned photometric and kinematic components, except 
for NGC 474 which exhibits obvious irregular shells. This contrasts 
with slow rotators which have significant kinematic misalignments 
and/or velocity twists. 

(vi) All slow rotators, besides the three slowest rotators (for 
which stellar velocities are consistent with zero everywhere), ex- 
hibit a kiloparcsec-scale stellar KDC (see also Paper VIII). 

These results are confirmed when including 18 additional 
galaxies observed with SAURON. Although our sample of 48 early- 
type galaxies is not a good representation of the galaxy luminosity 
distribution in the nearby Universe, it is the first to provide such a 
constraint for numerical simulations of galaxy formation and evo- 
lution. It is clear that a significantly larger and unbiased sample is 
required to probe the true distribution of A_r in early-type galax- 
ies and the fraction of slow and fast rotators. Even so, our sample 
should serve as a reference point for detailed numerical simulations 
of galactic systems in a cosmological context. 



6 CONCLUSIONS 

Using two-dimensional stellar kinematics, we have argued that 
early-type galaxies can be divided into two broad dynamical 
classes. We have devised a new parameter A_r to quantify the spe- 
cific angular momentum in the stellar component of galaxies. This 
parameter, as applied within one effective radius to the 48 E and SO 
galaxies of our SAURON sample, allowed us to disentangle the two 
classes of fast and slow rotators with a threshold value of 0.1. As 
emphasised in Sect. 13,21 fast and slow rotators exhibit qualitatively 
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APPENDIX A: THE ANGULAR MOMENTUM OF 
GALAXIES AND X R 

The criterion we wish to define should clearly separate the slow 
and fast rotators. The first guess is to take something equivalent to 
the total angular momentum, which should be expressed as some- 
thing like (J*l|V|) where X is the distance to the spin axis. To 
avoid having to determine that axis, we chose (7i|V|), where R 
is the distance to the centre. The absolute value of V is used as we 
are interested in the presence of local streaming motion (e.g. see 
NGC 4550). This is then naturally normalised by (Xy/V 2 + a 2 ) 
where V 2 + a 2 is the second order velocity moment. 

Another route would be to get closer to some physical quan- 
tity such as the angular momentum per unit mass. This can be rep- 
resented via the spin parameter (e.g. ?): 



A = 



GM 2 - 5 



(Al) 



where J, E, and M are the total angular momentum, energy and 
mass, and G is Newton's constant of gravity. The derivation of A 
in principle requires an accurate distance for each galaxy. Since we 
wish to stay close to our SAURON measurements we need to define 
something which can be readily evaluated. This is possible by first 
rewriting A as: 



A = 



J/Ms/\E\/M 
GM 



(A2) 



For two-dimensional integral-field data, we can approximate the 
spin parameter by using the scalar virial relations and sky aver- 
aging over the field-of-view, weighted with the surface brightness. 
The scalar virial relations are 2E — —2K — W, with the to- 
tal kinetic and potential energy related to the mean-square speed 
of the system's stars Vrms respectively as 2K/M = V£ MS and 
W/M — —GM/r g . Here, r g is the gravitational radius, which 
can be related to the half-mass radius rn given the mass model of 
the stellar system. For e.g. a spherical symmetric ? model, we have 
r g = 2r h . 

When relating r g (or ru) to the observed effective radius R e , 



we have to take into account the projection as a function of the 
viewing angle of the stellar system which in general is not spher- 
ical. Furthermore, while r g and are related to the mass dis- 
tribution, R e depends on the light distribution. Therefore, given 
r g = kr Re the conversion factor kr is a function of the mass 
model, viewing direction and mass-to-light ratio. In a similar way 
a conversion factor is needed in the approximation of the total an- 
gular momentum by the observed radius times velocity: J/M = 
kj R\ V| (see ?). Furthermore, while V^ MS is the total mean-square 
speed, we only observe velocity moments projected along the line- 
of-sight, e.g. V and a. In order to retrieve Vrms> we should there- 
fore allow multiplying factors kv and ks associated with the corre- 
sponding observed V 2 and a 2 (in the case of an isothermal sphere, 
kv = 1 and Ks = 3). 

We thus find the following expressions: 



J/M = kj (R\V\), 

2E/M = - (k v V 2 + kso- 2 ), 

GM = k r {R(k v V 2 + ks a 2 )) , 

resulting in an approximate spin parameter 

kj (R\V\) V(kvV 2 +k s a 2 } 



A 



krV2 (R{K V V 2 + Ksa 2 )) 



(A3) 
(A4) 
(A5) 

(A6) 



For simplicity, we will use the observed second order velocity mo- 
ment V 2 + a 2 , so that both kv — 1 and ks = 1, and thus define 
X f as: 



A 



(R\V\) ^(V 2 +a 2 ) 



(R(V 2 + a 2 )) ■ 
We finally approximate A f with A_r by writing 

(R\v\) 



{RVV 2 + a 2 } 



(A7) 



(A8) 



When deriving both the values of A / and Xr for the 48 E and SO 
SAURON galaxies, we indeed find a tight relation: 



Xr = (0.95 ± 0.04) x X f . 



(A9) 



This relation is in fact also valid for the two-integral models as de- 
rived in Appendix [B] so that we expect An to be roughly equal 
to (though slightly smaller than) Xf. For typical values of the 
conversion factors, kj = 2, kr — 3, we therefore find that 

A ~ a/2/3 X r . 

An alternative to Xr is the unbounded A^ = (R\V\)/{Ra). 
In the context of early-type galaxies we favour A_r which includes a 
mass-like normalisation by the second order velocity moment V 2 + 
a 2 . For isotropic (two-integral) models, this implies Xr ~ %/e. For 
slow rotators Xr ~ X' r, since these galaxies have by definition 
small mean stellar velocities. 



APPENDIX B: TWO-INTEGRAL MODELS 

In order to examine the behaviour of Xr in more detail, we con- 
structed two-integral dynamical Jeans models for 7 SAURON galax- 
ies using the MGE formalism, namely NGC 524, 3377, 4459, 4486, 
4552, 4621, and 5813. Details on the assumptions and resulting 
mass models for these galaxies can be found in Paper IV. For each 
galaxy, we assumed 8 different inclinations (5, 15, 25, 35, 50, 65, 
80, and 90°, the latter corresponding to an edge-on view), and de- 
rived the first two line-of-sight velocity moments up to 1 R e . 

These models were used to obtain measurements of Xr, Xf 
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and V/a, to determine whether these three quantities behave sim- 
ilarly with respect to basic parameters such as the inclination 
or anisotropy parameter. We first obtain that \r — (0.93 ± 
0.04) x \f, a correlation very similar to what was found for the 
observed 48 SAURON galaxies (see Appendix [A}- The variation 
of \r with inclination follows that of V/a within better than 
10% for two-integral models. These results do not depend on the 
global anisotropy, as long as the internal dynamics of the model 
is fixed before changing the viewing angle. A radial variation in 
the anisotropy profile can, however, significantly change this situ- 
ation, as \r includes an additional explicit dependence on radius. 
Assuming isotropy for all 7 MGE Jeans models, we observe a tight 
correlation between Xr and V/a. A simple approximation for the 
relation between Xr and V/a can then be obtained by introducing 
a scaling factor k: 



Xr = 



(RV) 



K (V/a) 



(RVV 2 + a 2 ) 



(Bl) 



(V/a) 2 



Using our two-integral models we find a best fit relation with k = 
1.2 ± 0.1, slightly higher than but still consistent with the best fit 
value hi = 1.1 i 0.1 obtained using our 48 SAURON Es and SOs. 
We can in fact directly solve for k in Eq. dBlt in terms of the values 
of Xr and V/a measured for individual two-integral models: there 
is in fact a rather large spread in the values of k, which ranges from 
1 and 1.4, depending on the galaxy. This shows that non-homology 
is an important driver of changes in Xr at constant V/a. 

Another important issue when deriving Xr and V/a comes 
from the presence of noise in the measurements. When V is suf- 
ficiently large, the random errors in the kinematic measurements 
will affect both Xr and V/a randomly. As V becomes small in 
absolute value, (V^ 2 ) and {R \V\) become increasingly biased (be- 
ing artificially increased), which affects estimates of Xr and V/a. 
In Fig. IB II we show the effect of a random error of 10 km s _1 in 
both V and a (typical in the outer part of SAURON kinematic maps; 
see Paper III) on both Xr and V/a, in terms of Xr. A is the dif- 
ference between the measured value (with noise) and the expected 
one (noiseless) and is plotted with respect to the true Xr value. The 
effect is in general relatively small, of the order of 7% on Xr and 
almost 15% on V/a, for Xr ~ 0.1. When V is close to zero every- 
where, as is the case for the three SAURON galaxies with the lowest 
Xr values, the presence of noise yields an overestimate of about 
0.02 in Xr (depending on the velocity dispersion values), an effect 
consistent with the observed values for these galaxies (see Fig.[5]l. 
As expected, Fig. lBll also shows that Xr is in general about a factor 
of two less sensitive to errors in the kinematic measurements than 
V/a. 
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Figure Bl. Effect of noise in the kinematics on the estimate of A^j (filled 
symbols) and V/cr (empty symbols). A, the difference between the mea- 
sured value (with noise) and the expected one (noiseless), is plotted against 
Xr. The input noise has been set to 10 km s _1 for both V and a. 



